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ABSTRACT 


The  Naval  Postgraduate  School  has  an  ongoing  effort  to  study 
geomagnetic  noise  and  micropulsations  in  the  ULF  frequency  range 
(.05  4  f  <  10  Hz).  Data  is  collected  by  three  orthogonally 
mounted  coils  at  a  remote  land  site  and  telemetered  to  the 
laboratory  for  computer  analysis. 

To  isolate  data  containing  a  micropulsation  event,  time 
series  plots  of  the  magnetic  field  were  generated.  The 
development  of  a  double  running  average  routine  made  possible 
the  isolation  of  micropulsations  in  large  data  sets.  A  type 
PC-3  micropulsation  was  found  and  the  coherence,  ellipticity  and 
polarization  properties  were  determined  as  follows:  Coherence  * 
0.99,  Degree  of  Polarization  =  0.99  and  the  micropulsation  was 
elliptically  polarized. 

Power  spectral  density  (PSD)  plots  summarizing  data  segments 
about  two  hours  long  occasionally  contained  structures  found  to 
be  artificial  and  not  representative  of  natural  phenomena  in  the 
geomagnetic  field.  Methods  to  avoid  this  anomalous  behavior  in 
PSD  plots  are  suggested.  '*r;? 
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I.  INTRODUCTION 


This  thesis  is  part  of  an  ongoing  effort  at  the  Naval 
Postgraduate  School  to  analyze  ULF  geomagnetic  noise  and 
micropulsations.  The  results  of  the  Naval  Postgraduate  School 
studies  could  impact  communications  systems  or  systems  in  which 
geomagnetic  noise  and/or  micropulsations  introduce  operational 
difficulties . 

An  objective  of  this  thesis  was  to  develop  the  software 
necessary  for  obtaining  computer  generated  plots  of  geomagnetic 
field  versus  time.  Then,  with  the  plots,  sections  of  data 
containing  obvious  geomagnetic  micropulsations  were  located. 
Those  sections  were  analyzed  to  determine  the  power  spectral 
density  (PSD),  coherence,  degree  of  polarization  and 
ellipticity  of  the  micropulsation  using  already  developed 
software  [Ref.  1]. 

After  analyzing  some  of  the  plots  (PSD,  coherence,  etc.) 
produced  by  the  software  documented  in  [Ref.  1],  it  was 
demonstrated  that  problems  existed  in  the  software.  New 
software  was  developed  to  obtain  times  series  raw  voltage  and 
geomagnetic  field  plots.  Validation  of  the  software  required 
an  on-site  experiment.  Documented  within  this  thesis  is  the 
software  to  produce  time  series  voltage  and  geomagnetic  field 
plots  and  the  questions  the  validation  experiment  raised 
concerning  the  sensing  system  and  software. 
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Using  the  newly  developed  software,  a  micropulsation  event 
was  located.  A  double  running  average  software  routine 
accentuated  the  micropulsation  to  permit  the  determination  of 
its  frequency  and  amplitude.  The  micropulsation  event  was 
further  analyzed  to  obtain  information  on  coherence,  degree  of 
polarization  and  ellipticity. 


II .  BACKGROUND 


A.  COLLECTION  SYSTEM 

The  Naval  Postgraduate  School  has  two  operating  geomagnetic 
sensing  systems.  The  land  site  is  located  at  the  La  Mesa 
Village  Housing  Area  and  has  three  orthogonally  mounted 
induction  coils.  The  second  operating  system  contains  two 
orthogonally  mounted  induction  coils  that  can  be  placed  on  the 
bottom  of  Monterey  Bay.  When  both  systems  are  operating,  the 
data  is  collected  simultaneously.  The  ground  work  has  been  laid 
for  a  third  system  located  at  Chew's  Ridge.  Three  orthogonal 
and  simultaneous  measurements  can  be  taken  only  between  two 
sites  since  the  School  has  only  five  induction  coils  and  a  sixth 
just  calibrated.  This  thesis  used  data  from  the  La  Mesa  Village 
site  and  the  software  associated  with  it. 

An  induction  coil  senses  the  geomagnetic  field  fluctuations 
based  on  Faraday's  Law  of  induction,  induced  emf  =S"N^‘»  which 
describes  the  relation  between  the  induced  emf  and  the  time 
rate  of  change  of  magnetic  flux  through  the  coil.  N  is  the 
number  of  turns  in  the  coil. 

The  voltage  induced  in  the  coils  is  sampled  32  or  64 
times/second  and  is  amplified  approximately  one  million  times. 
The  voltage  data  undergoes  pulse-code-modulation  (PCM)  for 
noiseless  VHP  radio  link  transmission  from  the  La  Mesa 
Village  site,  to  the  recording  system  in  the  Geophysics  Signal 
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Processing  Laboratory  located  in  Spanagel  531.  Data  reduction 
on  the  Naval  Postgraduate  School's  IBM  3033  computer  follows 
digitization  of  the  PCM  data.  A  block  diagram  of  the  system  is 
provided  in  Figure  2.1.  A  more  detailed  description  of  the 
system  is  given  in  Reference  2. 

1.  System  Calibration 

For  purposes  of  documentation,  the  calibration  of  the 
third  coil  for  either  the  Chew's  Ridge  or  Monterey  Bay  system 
is  being  included  in  this  thesis.  The  general  method  of 
induction  coil  calibration  will  be  covered  here.  For  a  more 
indepth  explanation,  see  Reference  1. 

The  sensing  coil  and  its  associated  electronics  are 
calibrated  using  a  Helmholtz  coil  apparatus  to  establish  a 
uniform  magnetic  field.  By  placing  the  sensing  coil  into  the 
field  produced  by  the  Helmholtz  coil,  one  can  establish  a 
relationship  between  the  Helmholtz  coil  magnetic  field  magnitude 
and  frequency  and  the  coil  system  voltage.  Once  the  response  of 
the  coil  system  is  known  for  the  frequency  band  of  interest,  it 
can  be  included  in  the  system  transfer  function  used  in  data 
reduction. 

The  experiment  to  determine  the  sixth  coil  system 
transfer  function  was  done  at  the  La  Mesa  Village  site.  A 
Wavetek  Model  142  signal  generator  supplied-  a  sinusoidal  current 
to  a  1.22  meter  diameter,  0.61  meter  high  Helmholtz  coil.  The 
current  was  measured  across  a  996  ohm  resistor  in  series  with 
the  Helmholtz  coil.  A  Hewlett  Packard  HP-3582A  spectrum 
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Figure  2.1.  System  Data  Flow 


analyzer  measured  the  output  voltage  of  the  coil  system  as  a 
function  of  frequency. 

The  experiment  measured  the  response  of  the  coil  system 
in  the  frequency  band  .05  hertz  to  20  hertz  for  applied  fields 
of  0.02,  0.2  and  1.0  nanoteslas.  The  coil  system  response  is 
shown  in  Figure  2.2.  The  transfer  function  of  this  coil  system 
corresponds  closely  to  those  of  coil  systems  previously 
calibrated.  The  only  noticable  difference  is  that  the  0.02 
nanotesla  response  at  15  hertz  is  approximately  15  volts/ 
nanotesla  higher  than  those  measured  before  [Ref.  1].  It 
should  be  noted  that  this  difference  will  not  influence  data 
analysis  in  the  frequency  range  of  interest.  The  transfer 
function  algorithm  is  listed  in  Table  2.1. 

B .  PREVIOUS  SOFTWARE 

Software  has  been  developed  (see  References  1  and  3),  to 
extract  power  spectral  densities  for  each  coil,  coherence 
between  the  coils,  degree  of  polarization  in  the  three  measure¬ 
ment  planes,  ellipticity  in  the  three  measurement  planes  and  a 
variety  of  Stokes  parameters.  The  following  is  a  basic 
description  of  the  computer  code  that  produces  the  above 
mentioned  plots.  First,  the  voltage  data  is  read  off  the 
digital  tape  using  a  subroutine  called  RD  provided  by  Or.  Tim 
Stanton  of  the  Oceanography  Department.  A  parameter  ISEC, 
representing  the  number  of  seconds  one  wishes  to  advance  the 
tape,  is  frequently  used.  ISEC  establishes  the  number  of 
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Table  2.1.  System  Transfer  Function  Algorithm 


DC  9  L  =  1  ,N 
FRQ*FR£Q IL) 

IF  ( FRQ .LE.25.IGO  TO  1 
XX(L)«XX(tT/28. 

GC  TC  3 

IF(FRQ.LE.15.JGG  TO  2 
XX(L)*XX<L1/(1Q5.5-3.14*FR3) 
YY(L)=YY(L)/  <m.32-7.5aa*FRQJ 
ZZ(L)*ZZ(L)/(177.26-7.434*FRQ) 

GC  TO  8 

IF<FRQ<LE< 10*) GO  TO  3 
XX(LJ*XX(Li/  (  5,9  58^  FRQ— 30.97 ) 
YY(L)*YY(L1/ (7.1 66*  FRQ-39 .99 i 
ZZ(U*ZZ(L)/(6.49*FRQ-32.35) 

GC  TO  3 

IF(FRQ.LE.7.5JGO  TO  4 
XX(L)=XX(L)/ (3.492*  FRQ-6. 31 ) 

YY  ( L)*YY  (L  )  /  (4.252*FRQ-i0.85i 
ZZ(L)*ZZ(Li/  (4.04**  FRQ-7  .89  > 

*jL  TC  3 

IF  <F RQ  •  L E.  5  •  J  Gu  TC  5 
XX(Li«XX(LJ/ <2.6311 *FRQ+0.i4oo 7 i 
YY(L)=YY(L)/(3.0  12* FRQ- 1 . 55 J 
ZZ( L)-ZZ (L) / (3.134* FRQ-i.44) 

GC  TO  8 

IF(FRQ.LE.3.  JGO  TO  6  ,  n  ,  „  „ 

XX(LJ«XX(L)/  <2.63ll*FRQ*0.U6*7) 
YY(L)*YY(U/  (2.702*  FRQ) 

ZUU-LZ  (LJ/  <Z.92*FRQJ 
GC  TO  8 

XX(LJ«XX(U/  <2.72*FRQI 
GC  TC  7 

continue 

TF(Li«(XX(L)*COSO  +  YY  (U  *CuSDl  )*COSoO 
CCNTINJE 


ZZ(L)*CGS30 
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times  RD  is  called  and  the  data  not  stored.  Previously,  tape 
advances  of  approximately  30  seconds  or  ISEC  =  30  were  -used. 
Once  the  tape  is  advanced  the  desired  amount,  data  analysis  can 
begin.  Each  digital  tape  contains  about  90  minutes  or  1.728  x 
10*  pieces  of  data  per  axis  resulting  from  the  sampling  rate  of 
32  samples/second.  Because  of  the  large  amount  of  data,  the 
analysis  has  to  be  accomplished  in  blocks.  Data  analysis  takes 
place  on  256  seconds  or  8192  frames  of  data  at  a  time.  RD  is 
called  8192  times  and  places  the  voltages  for  each  coil  in  an 
array.  The  voltages  are  integer  values  from  0  to  4096 
representing  -5  volts  to  +5  volts.  Ideally  the  value  2048 
represents  0  volts.  The  integers  between  0-4096  are  normalized 
to  +  5  volts  and  placed  in  new  complex  arrays.  These  arrays 
containing  8192  pieces  of  data  are  Fourier  transformed  to  the 
frequency  domain  where  the  system  transfer  function  is  applied. 
After  applying  the  transfer  function,  the  calculations  begin 
for  a  determination  of  PSD's,  coherence,  ellipticity,  etc.  for 
the  first  block.  After  they  are  complete,  the  next  block  of 
8192  pieces  of  data  starts  through  the  program.  Once  the 
desired  amount  of  data  has  been  processed,  a  NONIMSL  subroutine 
called  DRAWP  plots  the  data.  For  a  more  complete  description 
of  this  software,  see  Reference  1,  Reference  3  and 
Reference  4. 
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C.  ANALYSIS  OP  INITIAL  DATA  PLOTS 

Many  of  the  PSD's,  coherences,  ellipticities  and  degree  of 
polarization  plots  produced  by  the  software  of  Reference  1  and 
Reference  3  have  unexpected  characteristics.  The  PSD's  show 
“humps"  that  are  shown  in  Figures  2.3-2. 5.  These  "humps"  are 
not  characteristic  in  "normal"  PSD's.  The  coherence  plots  show 
a  coherence  of  essentially  one  degrading  to  hash  in  the 
frequencies  greater  than  2  hertz.  Figures  2. 6-2. 8.  We  see  a 
similar  behavior  in  the  plots  of  degree  of  polarization.  These 
plots  indicate  a  very  highly  polarized  field  from  0.02  hertz  to 
1.0  hertz.  Figures  2.9-2.11.  The  ellipticity  plots  show  a  very 
linearly  polarized  field  for  these  frequencies  in  Figures 
2.12-2.14. 

After  examining  the  time  series  magnetic  field  plots,  we 
noted  a  section  of  very  large  magnitude  fluctuations  (about  one 
to  three  orders  of  magnitude  greater  than  the  general  back¬ 
ground)  at  the  beginning  of  the  digital  tapes  Figures  2.15-2.17. 
We  found  that  if  the  digital  tape  was  advanced  past  these  large 
fluctuations  the  output  (PSD,  coherence,  etc.)  appeared 
"normal".  Figures  2.18-2.29.  If  the  program  was  executed  on 
just  the  section  containing  the  large  fluctuations,  the  PSD's 
had  the  "humped"  characteristic,  the  coherences  were  one,  the 
degree  of  polarization  was  very  high  and  the  ellipticity  was 
linear.  Figures  2.30-2.41.  This  demonstrates  that  these  large 
fluctations  dominate  the  output  of  an  analysis  in  the  frequency 
domain.  A  tape  advance  of  300  seconds  or  ISEC  ■  300  is 
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(Hz)  (0.5  units/in),  16  Averages 
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Figure  2.8.  Coherence  Z  and  Z  Coils,  17  August  82,  2240-2348  Local. 

Coherence  Z  and  X  Coils  (0.2  units/in)  vs.  Log  Frequency 
(Hz)  (0.5  units/in),  16  Averages 


Figure  2.9.  Degree  of  Polarization  X-Y  Plane,  17  August  82,  2240-2348  Local 
Degree  of  Polarization  (0.2  units/in)  vs.  Log  Frequency  (Hz) 
(0.5  units/in),  16  Averages 


equency 


Figure  2.11.  Degree  of  Polarization  Z-X  Plane,  17  August  82,  2240-2348  Local 
Degree  of  Polarization  (0*2  units/in)  vs.  Log  Frequency 
(Hz)  (0.5  units/in)  ,  16  Averages 


Figure  2.12.  Ellipticity  X-Y  Plane,  17  August  82,  2240-2348  Local. 

Ellipticity  (0.5  unit/in)  vs.  Log  Frequency  (Hz)  (0.5  unit/in) 


Figure  2.14.  Ellipticity  Z-X  Plane,  17  August  82,  2240-2348  Local. 

Ellipticity  (0.5  units/in)  vs.  Log  Frequency  (Hz)  (0.5  units/in) 


Figure  2.15.  X-Coil  Magnetic  Field,  17  August  82,  2240-2348-  Local. 

Amplitude  (nanoteslas  :  10  units/in)  vs.  Time  (seconds  :  500  units/in) 


Figure  2.16.  Y-Coil  Magnetic  Field,  17  August  82,  2240-2348  Local. 

Amplitude  (nanoteslas  :  10  units/in)  vs.  Time  (seconds  :  500  units/in) 


Figure  2.17.  Z-Coil  Magnetic  Field,  17  August  82,  2240-2348  Local. 

Amplitude  (nanoteslas  :  5  units/in)  vs.  Time  (seconds  :  500  units/in) 
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Coils  (0.2  units/in)  vs.  Log  Frequency  (Hz) 
,  16  Averages 


Figure  2.23.  Coherence  Z-X  Coils,  17  August  82,  2245-2353  Local. 

Coherence  Z-X  Coils  (0.2  units/in)  vs.  Log  Frequency  (Hz) 
(0.5  units/in),  16  Averages 


Figure  2.26.  Degree  of  Polarization  Z-X  Plane,  17  August  82,  2245-2353  Local 
Degree  of  Polarization  (0.2  units/in)  vs.  Log  Frequency  (Hz) 
(0.5  units/in),  16  Averages 


Figure  2.27.  Ellipticity  X-Y  Plane,  17  August  82,  2245-2353  Local 
Ellipticity  (0.5  units/in)  vs.  Log  Frequency  (Hz) 
(0.5  units/in),  16  Averages 


Figure  2.28.  Ellipticity  Y-Z  Plane,  17  August  82,  2245-2353  Local. 

Ellipticity  (0.5  units/in)  vs.  Log  Frequency  (Hz)  (0.5  units/in) 
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Degree  or  Polarization  X-Y  Plane,  17  August  82,  2240-2250  Loca 
Degree  of  Polarization  (0.2  units/in)  vs.  Log  Frequency  (IIz) 
(0.5  units/in),  2  Averages 
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Figure  2.38.  Degree  of  Polarization  Z-X  Plane,  17  August  82,  2240-2250  Local 
Degree  of  Polarization  (0.2  units/in)  vs.  Log  Frequency  (Hz) 
(0.5  units/in),  2  Averages 


Figure  2.39.  Ellipticity  X-Y  Plane,  17  August  82,  2240-2250  Local 
Elliptic! ty  (0.5  units/in)  vs.  Log  Frequency  (Hz)  (0 


Figure  2.41.  Ellipticity  Z-X  Plane,  17  August  82,  2240-2250  Local 
Ellipticity  (0.5  units/in)  vs.  Log  Frequency  (Hz) 
(0.5  units/in),  2  Averages 


recommended  at  the  beginning  of  the  programs  producing  the  above 
mentioned  plots.  Time  series  plots  should  also  be  made  to 
provide  a  further  check. 

An  explanation  of  the  source  of  the  large  fluctuations  is 
that  the  PCM  to  digital  conversion  begins  before  the  actual  PCM 
data  starts  on  the  tape.  The  decoding  equipment  considers  the 
noise  on  the  section  of  tape  before  the  PCM  data  starts  as 
signal. 

The  above  analysis  clearly  proves  that  the  "hump"  like 
structure  in  the  PSD's  and  the  “nice"  behavior  observed  in  the 
corresponding  coherence,  ellipticity  and  degree  of  polarization 
are  artifical  and  do  not  represent  any  real  phenomena  in  the 
geomagnetic  field.  Thus,  the  data  in  Figures  2.3-2.14  and 
Figures  2.30-2.41  does  not  appear  to  be  reliable.  Unfortu¬ 
nately,  this  artificially  induced,  and  therefore  unphysical, 
behavior  of  the  polarization  has  been  published  recently 
[Ref.  7].  The  conclusion  found  in  Reference  7  that 

Geomagnetic  fluctuations  in  land  are  well  polarized 
below  1.0  Hz  and  that  the  polarization  is  quite 
frequently,  though  not  consistently,  linear.  The 
situation  under  the  sea  is  similar.  Furthermore, 
a  land-sea  coherence  study  in  the  horizontal  plane 
showed  very  high  coherence  below  1.0  Hz. 

are  not  substantiated  by  the  physical  data.  Plots  of  the 

artificially  induced  behavior  are  published  in  References  1,  2, 

and  3. 
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The  4  Hz  peak  present  in  the  PSD's  shown  in  Figures  2.18- 
2.20  is  the  result  of  60  Hz  aliasing.  60  Hz  magnetic  fields 
generated  by  local  power  lines  exist  at  the  measurement  site. 
Digitization  at  a  32  Hz  or  64  Hz  rate  of  the  analog  voltages 
from  the  sensing  coils  shifts  the  true  frequency  of  the  power 
line  signals  down  at  4  Hz.  The  60  Hz  aliasing  results  in  the 
anomalis  behavior  at  4  Hz  seen  in  the  coherence,  degree  of 
polarization  and  ellipticity  plots  shown  in  Figures  2.21-2.29. 

Another  peak  is  present  at  about  1.2  Hz  for  the  PSD's  in 
Figures  2.18-2.20.  For  this  data  set,  the  1.2  Hz  signal 
correlates  with  a  1.2  Hz  signal  identified  in  the  time  series 
data.  The  origin  of  the  1.2  Hz  signal  is  unknown.  The  4  Hz  and 
1.2  Hz  peaks  are  of  little  concern  in  this  study  since  both  are 
well  above  the  frequency  of  the  micropulsation  analyzed. 

Another  "suspicious"  feature  of  the  PSD's  is  the  large 
amount  of  hash  appearing  from  2  to  10  hertz.  We  will  gain 
further  insight  into  this  problem  in  analyzing  the  time  series 
voltage  data. 

From  the  analysis  of  the  abnormal  behavior  in  the  PSD’s, 
coherence,  etc.,  future  coil  data  sets  should  receive  a 
balanced  treatment  between  the  time  domain  and  the  frequency 
domain  to  avoid  misinterpretation. 


III.  TIME  SERIES  DATA 


A.  TIME  SERIES  VOLTAGE  DATA 
1.  Voltage  Software 

The  importance  of  time  series  data  as  a  check  became 
evident  when  we  uncovered  the  problem  with  the  initial  PSD's, 
coherences,  etc.  Computer  code  was  developed  to  produce  both 
unsmoothed  and  smoothed  times  series  voltage  plots.  The  reason 
for  producing  both  smoothed  and  unsmoothed  voltage  plots  will 
become  apparent  in  the  time  series  voltage  analysis  section. 
Time  series  voltage  plots  could  provide  a  means  of  monitoring 
the  data  for  system  "glitches"  which  might  later  produce 
erroneous  PSD's,  etc.,  and  suspicous  features  on  time  series 
magnetic  field  plots. 

The  Fortran  computer  program  that  produces  the  voltage 
plots  is  compiled  and  run  on  the  Naval  Postgraduate  School  IBM 
3033  VM  or  "batch"  system.  Both  the  unsmoothed  voltage  program 
(VOLTR )  and  the  smoothed  voltage  program  (VOLTS)  require  2.048 
megabytes  of  core  memory  and  approximately  15  minutes  of 
central  processor  unit  (CPU)  time  to  run.  These  figures 
provide  a  comfortable  margin  to  insure  that  the  program  does 
not  "bomb"  due  to  lack  of  core  or  run  time. 

The  design  of  both  VOLTR  and  VOLTS  is  the  same  with  the 
exception  that  a  smoothing  algorithm  is  applied  to  the  data  in 
the  case  of  VOLTS.  Both  programs  have  incorporated  in  them  the 
digital  tape  advance  package  using  ISEC  as  the  number  of 
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seconds  one  wishes  to  advance  the  tape.  Data  is  read  off  the 
digital  tape  in  blocks  of  8192  frames  by  the  subroutine  RD 
referenced  earlier.  As  before,  this  data  is  in  integer  form 
and  represents  voltages  between  ±5  volts  by  numbers  between  0 
and  4096.  The  integer  representing  0  volts  is  ideally  2048, 
however,  I  found  that  the  "real"  zero  values  wander.  The 
reason  for  this  "wander"  is  not  understood.  For  good  results, 
the  following  values  worked  well  for  the  respective  coils. 

X-coil  :  zero  «  1966 

Y-coil  :  zero  =  2085 

Z-coil  :  zero  =  2539 

The  voltage  values  normalized  to  +5  volts  are  placed  in  new 
real  arrays .  At  this  point  the  next  block  of  data  is  read  off 
the  digital  tape  and  goes  through  the  same  process.  VOLTR  and 
VOLTS  analyze  eight  blocks  or  34  minutes  of  data.  The  linear 
arrays  holding  the  voltage  data  are  dimensioned  to  65,536 
elements  given  a  sampling  rate  of  32  samples/ second.  For  every 
one  of  the  65,536  voltage  values  there  is  a  time  value  in 
another  array.  Thus  four  arrays  are  required  to  handle  the 
data,  three  arrays  containing  voltage  data  and  one  that 
contains  the  respective  times  for  each  of  the  measurements.  At 
this  point,  the  VOLTR  program  calls  the  NONIMSL  subroutine 
DRAWP  to  plot  the  data.  The  VOLTS  program  uses  a  smoothing 
algorithm  to  reduce  the  amount  of  high  frequency  hash. 

The  smoothing  algorithm  is  a  double  running  point 
average  executed  on  all  three  axes  of  data.  Since  1  was 
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interested  in  periods  of  approximately  10-45  seconds,  I 
averaged  over  144  points.  With  a  running  average  over. 23 
percent  of  the  oscillation,  the  unwanted  background  can  be 
removed  without  destroying  the  oscillations.  Averages  greater 
than  50  percent  of  the  period  would  partially  remove  the 
oscillation  of  interest.  The  construction  of  the  algorithm 
makes  changing  the  number  of  points  averaged  over  trivial. 

One  need  only  change  a  loop  index  and  three  divisors  to  obtain 
a  different  average.  The  smoothing  algorithm  is  provided  as 
Table  3.1.  The  computer  code  for  both  VOLTR  and  VOLTS  is 
attached  in  Appendix  A  and  Appendix  B. 

2.  Time  Series  Voltage  Data  Analysis 

One  can  gain  a  great  deal  of  insight  into  the  contents 
of  the  data  by  looking  at  the  smoothed  and  unsmoothed  voltage 
products.  Figures  3. 1-3. 3  are  unsmoothed  voltage  plots 
representing  the  first  34  minutes  of  the  digital  data  tape  GMDT 
11  recorded  17  August  82,  1301-1408  local  time.  A  one  minute 
tape  advance  was  used  in  producing  this  data.  Immediately  one 
notices  the  large  amount  of  60  hertz  aliasing  in  the  time  series 
unsmoothed  voltage.  The  60  hertz  aliasing  appears  as  4  hertz 
because  of  the  32  sample/second  sampling  rate.  The  sensing 
system  has  a  low  pass  filter  that  should  be  improved  to  remove 
the  aliasing  problem.  The  4  hertz  background  can  be  seen  more 
clearly  in  Figures  3. 4-3. 6.  Above  the  60  hertz  aliasing 
background,  we  see  voltage  spikes  ranging  in  magnitude  from  1 
volt  to  over  4  volts.  Expanding  the  time  scale  to  50 
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Table  3.1.  Smoothing  Algorithm 


DC  73  L2*lr2 
S*Q 

OC  74  IS*1 *65318 

SUMX«0.0 

SUM Y*0 • 0 

SUMZ*0.0 

CO  75  J*  Li  144 

SUMX«2ZXiIq+J)+SUMX 

aaafJi  wwfflf 
75  mmi**™*. 

ifilHitttfflHft:*: 

i*o+i 

74  CONTINUE 
73  CONTINUE 
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Amplitude  (volts  :  2  units/in)  vs.  Time  (seconds  :  500  units/in) 


Figure  3.4.  X-Coil  Voltage,  17  August  82,  1322-1323  Local. 

Amplitude  (Volts:  2  units/in)  vs.  Time  (seconds  :  10  units/in) 


seconds  of  data,  the  structure  of  the  voltage  spikes  becomes 
apparent  in  Figures  3.4-3. 6.  The  "square  well"  shape  i-s 
characteristic  of  PCM  dropouts  or  instances  during  the  PCM  to 
digital  conversion  synchronization  is  lost  between  the  tape 
recorder  containing  the  PCM  tape  and  the  decoding  unit.  We  do 
not  know  whether  PCM  dropouts  always  occur  on  all  three 
channels,  but  current  evidence  suggests  dropouts  are 
simultanoues  on  all  channels. 

Figures  3.7-3. 9  again  show  the  first  34  minutes  of  GMDT 
11  only  smoothed  by  the  144  point  running  average.  The  smooth¬ 
ing  algorithm  removes  the  60  hertz  alising.  The  PCM  dropouts 
appear  as  "glitches"  in  the  otherwise  quasi  sinusoid  charac¬ 
teristic  of  a  micropulsation. 

We  should  note  that  the  unusually  large  amount  of  hash 
present  between  2  hertz  and  10  hertz  in  the  PSD's  and  the 
average  slope  of  the  PSD  curve  for  those  frequencies  may  in  part 
be  due  to  the  PCM  dropouts  that  occur  on  the  tapes .  A 
"deglitching"  algorithm  could  be  developed  to  remove  the  sudden 
large  amplitude  structures  injected  by  the  system  into  the 
voltage  data  if  hardware  modifications  fail. 

B.  TIME  SERIES  MAGNETIC  FIELD  DATA 
1.  Magnetic  Field  Software 

The  time  series  magnetic  field  software  is  basically 
the  same  as  that  which  produces  the  PSD’s,  etc.  [Ref.  1  and 
Ref.  3].  The  differences  begin  after  the  forward  Fourier 
transform  and  the  application  of  the  transfer  function.  We 
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Figure  3.7.  X  Coil  Voltage,  17  August  82,  1302-1336  Local. 

Amplitude  (volts  :  0.2  units/in)  vs.  Time  (seconds  :  500  units/in) 


Figure  3.9.  Z  Coil  Voltage,  17  August  82,  1302-1336  Local. 

Amplitude  (volts  :  0.1  units/in)  vs.  Time  (seconds  :500  units/in) 


iplitude  (nanoteslas  :  5  units/in)  vs.  Time  (seconds  :  500  units/in) 
point  smoothing 


point 


Amplitude  (nanoteslas  :  5  units/in)  vs.  Time  (seconds  :  500  units/in) 
144  point  smoothing 


Figure  3.12b.  Z  Coil  Magnetic  Field,  17  August  82,  1302-1336  Local, 

Amplitude  (nanoteslas  :  1  unit/in)  vs.  Time  (seconds  :  500  units/in) 
144  point  smoothing 


should  note  that  the  transfer  function  as  designed  ignores 
phase  information  since  it  is  just  a  least  squares  fit- to  data 
like  that  of  Figure  2.2.  Once  the  transfer  function  is 
applied,  the  reverse  transform  is  performed.  We  obtain  the 
complex  absolute  values  of  the  magnetic  field  data  (in 
nanoteslas)  and  place  the  data  in  new  real  arrays.  To  insure 
that  there  are  no  discontinuities  between  blocks  of  data. 

Figures  3.10a,  3.11a  and  3.12a,  all  the  blocks  are  connected  to 
the  first  by  adding  or  subtracting  a  constant  to  each  element, 
Figures  3.10b,  3.11b,  3.12b.  The  data  now  undergoes  the  same 
144  point  double  running  average  as  the  voltage  to  bring  out 
fluctuations  with  periods  of  10-45  seconds.  Again  DRAWP  is  used 
to  plot  the  data.  The  computer  program  LVFTC1  produces  the 
smoothed  magnetic  field  plots.  It  is  provided  in 
Appendix  C. 

2.  Smoothed  Magnetic  Field  Data  Analysis 

Again  looking  at  the  first  34  minutes  of  digital  data 
tape  GMDT  11,  we  see  in  the  magnetic  field  data  the  same 
exceptional  structures  as  seen  in  the  voltage  data,  Figures 
3.10-3.12.  The  glitches,  however,  do  highlight  one  problem 
rather  well.  From  Figures  3. 7-3. 9  we  note  that  positive 
voltage  spikes  produce  negative  magnetic  field  spikes.  This 
would  indicate  a  pi  phase  shift  is  introduced  either  in  the 
transfer  functions  or  in  the  digital  Fourier  transform.  After 
investigating  another  digital  tape  GMDT  3A,  18  August  82, 
0121-0251  local,  we  see  another  behavior.  Figures  3.13-3.15 
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Figure  3.13.  X  Coil  Voltage,  17  August  82,  2241-2315  Local. 

Amplitude  (volts  :  0.1  units/in)  vs.  Time  (seconds  :  500  units/in) 


Figure  3.14.  Y  Coil  Voltage,  17  Voltage  82,  2241-2315  Local. 

Amplitude  (volts  :  0.1  units/in)  vs.  Time  (seconds  :  500  units/in) 


are  the  smoothed  voltages  for  the  first  34  minutes  of  data 
on  GMDT3A  after  a  one  minute  tape  advance.  Figures  3.-16-3.18 
are  the  smoothed  magnetic  fields  for  the  same  time  period.  The 
X-coil  voltages  show  three  positive  spikes  while  the  X-coil 
magnetic  field  shows  three  negative  spikes .  This  supports  the 
theory  of  a  pi  phase  shift.  However,  the  Y-coil  voltage  shows 
two  positive  spikes  and  one  negative  spike,  while  the  Y-coil 
magnetic  field  has  three  positive  spikes.  Thus  it  appears  we 
have  a  rather  arbitrary  phase  problem  that  may  find  its  origin 
in  the  transfer  functions  that  do  not  maintain  the  correct 
phase . 

From  Figures  3.16-3.18  we  note  exceptional  features 
that  do  not  correspond  to  voltage  spikes.  The  exceptional 
features  are  separated  by  time  periods  of  256  seconds  which  is 
the  size  of  the  blocks  analyzed.  These  features  persist  with 
the  data  block  connection  scheme  and  144  point  smoothing.  Their 
origin  is  not  understood.  Some  data  blocks  look  like  Figure 
3.19  with  the  end  points  having  about  the  same  value.  However, 
some  have  end  points  as  illustrated  in  Figure  3.20.  The  varying 
end  point  values  are  also  shown  in  Figures  3.10a,  3.11a  and 
3.12a.  This  may  be  a  source  of  the  large  amount  of  hash  between 
2  hertz  and  10  hertz  seen  in  the  PSD's. 

After  looking  at  the  magnetic  field  data  it  was  obvious 
a  validation  experiment  was  needed  for  the  sensing  system  and 


software 


BLOCK  ENDS 


X  Coil  Magnetic  Field,  17  August  82,  2241-2315,  Local. 

Amplitude  (nanoteslas  ;  1  unit/in)  vs.  Time  (seconds  :  500  units/in) 


Figure  3.17.  Y  Coil  Magnetic  Field,  17  August  82,  2241-2315  Local. 

Amplitude  (nanoteslas  :  0.5  units/in)  vs.  Time  (seconds  :  500  units/in) 


PCM  DROPOUTS 


Figure  3.18.  Z  Coil  Magnetic  Field,  17  August  82,  2241-2315  Local. 

Amplitude  (nanoteslas  :  2  units/in)  vs.  Time  (seconds  :  500  units/in) 


IV.  SYSTEM  AND  SOFTWARE  VALIDATION 


A.  EXPERIMENTAL  APPARATUS 

The  inconsistencies  in  phase  between  the  time  series 
voltage  and  magnetic  field  plots  prompted  an  experiment  to 
validate  the  sensing  system  and  software.  The  sensing 
apparatus  consisted  of  the  X-coil  and  a  Schoenstadt  fluxgate 
magnetometer.  The  sensing  axis  of  both  sensors  were  horizontal 
and  oriented  North-South.  Figure  4.1  shows  the  test  set-up. 
Instead  of  an  on  axis  arrangement,  the  set-up  shown  was  used 
because  of  limited  resources.  The  fluxgate  magnetometer  is 
closer  to  the  source  since  it  has  less  sensitivity  at  the 
frequency  of  interest  than  the  coils.  Recorded  on  the  usual 
X-channel  was  the  output  from  the  X-coil.  The  Y-channel 
contained  data  from  the  coil  after  passing  through  an  additional 
Krohn-Hite  model  3321  lowpass  filter  with  a  10  hertz  cutoff.  On 
the  Z-channel  was  the  output  from  the  fluxgate  magnetometer. 

All  three  channels  were  recorded  simultaneously  and  a  64 
sample/second  sample  rate  was  used.  An  on-site  chart  recorder 
monitored  the  fluxgate  and  double  filtered  coil  voltage  output 
real  time.  The  experiment  produced  two  digital  data  tapes  each 
approximately  46  minutes  in  length  denoted  GMTT1A  and  GMTT1B. 

Artificial  magnetic  fields  were  introduced  into  the  region 
by  applying  sinusoidal  currents  from  a  Wavetek  signal  generator 
to  a  cylinder  (diameter  *  0.13m)  wrapped  with  150  turns  of 
wire.  Location  of  the  source  is  noted  in  Figure  4.1.  The 
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applied  sinusoidal  currents  varied  in  frequency  from  0.005 
hertz  to  10  hertz.  For  each  frequency,  two  magnitude  currents 
were  applied  to  the  source  to  demonstrate  a  linear  response  with 
amplitude. 

The  fluxgate  magnetometer  output  is  correct  in  phase  and 
amplitude  to  at  least  1  hertz.  The  transfer  function  between 
voltage  and  magnetic  field  for  the  fluxgate  magnetometer  is  a 
constant  of  10  nanoteslas/volt.  Thus  the  fluxgate  provides  an 
excellent  check  for  the  coil  transfer  functions,  the  RD  data 
stripping  subroutine  and  the  digital  fourier  transform. 

B .  EXPERIMENT  RESULTS 

1.  Voltage 

From  the  real  time  chart  record  of  voltage  for  the 
fluxgate  and  double  filtered  coil,  we  can  check  the  computer 
generated  voltages  using  VOLTR  and  VOLTS.  Table  4.1  shows  the 
chart  record  and  computer  generated  voltages  for  the  fluxgate 
and  double  filtered  coil  at  various  frequencies.  Only  the 
large  amplitude  oscillations  are  measured  for  the  double 
filtered  coil  because  of  poor  signal- to-noise,  for  the  small 
amplitude  oscillations. 

Table  4.1  shows  good  agreement  between  the  chart  record 
and  computer  output  for  the  double  filtered  coil.  However, 
the  fluxgate  magnetometer  voltages  undergo  a  constant 
transformation  of  0.831  from  chart  record  to  computer  plot. 

This  may  be  explained  by  a  problem  in  the  data  stripping 
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time  and  computer  generated. 


algorithm.  The  present  agreement  between  the  chart  record  and 
computer  generated  voltages  is  shown  in  Table  4.2. 

Another  more  serious  problem  became  obvious  during  the 
examination  of  the  computer  generated  voltage  plots.  For 
approximately  the  first  20  minutes  of  data  recording  the 
fluxgate  magnetometer  was  not  connected  to  the  PCM  board. 
Figures  4. 2-4. 4  show  the  two  respective  sensing  devices. 

Figures  4.2  and  4.3  both  contain  "real"  signal  while  Figure 
4.4  does  not  contain  "real"  signal  until  the  amplitude 
discontinuity  indicating  connection  of  the  fluxgate 
magnetometer  to  the  PCM  board.  Figures  4. 2-4.4  indicate  that 
"pick-up"  is  occurring  probably  between  the  channels  at  the  PCM 
board.  The  amplitude  of  the  cross-talk  voltage  is  about  0.15 
volts  peak-to-peak  at  the  frequency  of  0.1  hertz.  The 
amplitude  of  the  "real"  signal  on  the  double  filtered  coil  is 
0.29  volts  at  0.1  hertz.  Thus  51  percent  of  the  "real"  signal 
is  coupled  to  the  disconnected  channel.  Coupling  between 
channels  of  this  magnitude  raises  serious  questions  about  the 
computer  generated  data,  especially  for  coherence.  Discovering 
the  amount  of  the  channel  pick-up  if  the  channel  is  loaded 
requires  further  experimentation. 

Figures  4. 5-4. 7  show  an  interesting  feature  of  the  coil 
data  versus  fluxgate  magnetometer  data.  On  the  channels 
containing  the  coil  output  we  see  noise  on  top  of  the  signal. 
This  noise  is  not  present  on  the  channel  containing  the 
fluxgate  magnetometer  data.  Since  all  three  channels  passed 
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Table  4.2.  Percent  Agreement  Between  the  Chart  Record 
and  Computer  Generated  Voltages 
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through  the  PCM  encoder  and  decoder,  the  noise  is  probably  in 
the  coil  amplifiers  or  the  coils  themselves.  An  experiment  to 
investigate  the  true  origin  of  the  noise  should  be  done. 

2.  Magnetic  Field 

Given  the  dipole  equation  [Ref.  6]  and  the  experimental 
arrangement  in  Figure  4.1,  the  ratio  of  magnitude  of  the 
magnetic  field  seen  by  the  fluxgate  to  that  seen  by  the  coil 
can  be  calculated.  From  the  geometry  in  Figure  4.1  the  field 
at  the  coil  would  be  a  factor  of  8  less  than  that  seen  by  the 
fluxgate.  This  is  because  the  source  is  a  dipole  and  the  field 
of  a  dipole  falls  off  as  l/R  .  Detailed  consideration  of  the 
geometry  gives  a  ratio  of  approximately  8.09.  Table  4.3 
contains  the  magnetic  field  measured  at  the  fluxgate 
magnetometer.  Using  the  ratio  calculate^  from  the  geometry,  the 
expected  magnetic  field  amplitude  at  the  coil  is  listed  in 
Table  4.3. 

A  modified  form  of  computer  code  LFVTC1  generated  the 
magnetic  field  versus  time  plots  for  the  test  tapes.  Table  4.3 
shows  the  computer  generated  peak-to-peak  magnetic  field 
values  read  from  the  computer  plots.  From  Table  4.3  it  appears 
that  the  coil  magnetic  field  amplitude  is  in  error  by  a  factor 
of  approximately  three  (with  exception  to  the  low  frequencies). 
One  can  conclude  that  a  problem  exists  in  the  software. 

Figures  4. 8-4. 9  show  the  magnetic  field  data  for  the 
second  seventeen  minutes  of  digital  data  tape  GMTT1A.  Figure 
4.9  is  the  magnetic  field  measured  by  the  fluxgate  magnetometer 
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and  Figure  4.8  ia  the  magnetic  field  seen  by  the  coil.  For  the 
coil  we  see  discontinuities  every  256  seconds  corresponding  to 
the  block  size  of  data.  This  indicates  that  the  problems  in  the 
data  block  end  points,  as  mentioned  in  Section  III,  do  exist. 

It  is  not  clear  how  to  deal  with  this. 

The  first  300  seconds  of  Figures  4. 8-4. 9  is  expanded  in 
Figures  4.11-4.12.  With  the  expanded  data  a  phase  difference 
between  coil  and  fluxgate  can  be  determined.  By  checking  the 
phase  differece,  we  can  determine  whether  the  phase  independent 
transfer  function  is  introducing  error.  For  the  0.1  hertz  field 
plotted  in  Figures  4.10-4.11  the  coil  lags  the  fluxgate  by 
approximately  72  degrees.  In  Figures  4.12-4.13  we  see  that  at 
0.01  hertz  the  coil  sensitivity  has  fallen  and  a  reliable  phase 
difference  between  the  fluxgate  and  coil  cannot  be  made. 

Figure  4.14  shows  the  relationship  between  phase  angle  and 
frequency. 

The  above  experiment  shows  the  advantages  of  having  a 
real  time  "ground  truth"  record.  With  this,  one  can  avoid 
"working  in  the  blind".  Computer  products  will  have  a  check. 
Thus,  a  multi-channel  chart  recorder  should  be  purchased  to 
graph  in  real  time  the  x,  y,  and  z  inputs  to  the  PCM  board.  In 
addition,  electronic  circuitry  should  be  developed  and 
incorporated  so  that  the  chart  records  represent  at  least  a  low 
pass  filtering  of  the  coil  voltages.  Ideally,  conversion  of 
voltage  to  magnetic  field  could  occur  in  real  time  at  this  point 
given  the  appropriate  electronics. 
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V.  MICROPULSATIONS 


A.  THEORY 

Currently  there  is  no  theory  that  completely  models  the 
observed  micropulsations  and  geomagnetic  noise.  It  is 
generally  accepted  that  the  generating  mechanism  for 
micropulsations  is  an  interaction  between  the  solar  wind  and 
the  earth's  magnetosphere  [Ref.  5].  The  observed 
micropulsations  are  classified  by  their  period  and  regularity 
[Ref.  5].  Table  5.1  shows  the  micropulsation  classification. 

The  smoothing  algorithm  applied  to  the  magnetic  field  data  is 
designed  to  highlight  micropulsations  with  periods  between  10-45 
seconds.  Specifically  we  are  interested  in  micropulsations 
classified  as  Pc  3.  Pc  3  micropulsations  are  generally  regular 
and  have  a  more  or  less  distinct  period.  Thus  in  the  PSD  plots 
we  will  be  able  to  identify  the  micropulsation  by  an  apparent 
anomaly  at  its  frequency.  Micropulsations  have  amplitudes  of 
approximately  0.5  nT.  This  small  amplitude  is  very  small 
compared  to  a  main  magnetic  field  value  of  50,000  nT,  hence  the 
term  "micro"  pulsations.  Geomagnetic  fluctuations  are  a 
combination  of  micropulsations  and  a  general  background  noise. 
Micropulsation  events  grow  out  of  the  ever  present  background. 
The  time  series  magnetic  field  plots  will  provide  a  means  of 
obtaining  the  frequency  of  the  micropulsation. 
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B.  MICRO PULSATION  DATA  ANALYSIS 


The  first  1200  seconds  of  data  after  a  60  second  tape 
advance  on  digital  data  tape  GMDT  11,  17  August  82,  1301-1408 
Local,  may  contain  a  micropulsation  event.  Figures  3.10-3.12 
show  the  first  34  minutes  of  GMDT  11  after  a  60  second  advance. 
To  determine  whether  the  regular  pulsations  we  see  in  the  first 
third  of  Figures  3.10-3.12  (disregarding  the  suspicious  "cusps") 
is  a  micropulsation,  the  scale  must  be  expanded.  Figures  5.1- 
5.3  show  the  regular  pulsations  on  an  expanded  scale.  From 
Figures  5. 1-5. 3  the  period  of  the  oscillation  is  calculated  to 
be  17.2  seconds.  Thus  this  section  of  data  contains  a  type  Pc  3 
micropulsation. 

Executing  the  computer  code  LANDXYZ  [Ref.  1]  on  this 
section  of  data,  we  generate  the  PSD,  coherence,  degree  of 
polarization,  and  ellipticity.  Table  5.2  contains  the 
amplitude,  frequency,  coherence,  etc.  for  the  micropulsation. 

At  the  frequency  of  the  micropulsation  one  would  expect  a 
good  coherence,  -degree  of  polarization  and  ellipticity  since 
the  coil  system  is  oriented  arbitrarily  with  each  coil  sensing 
the  micropulsation.  The  observed  micropulsation  exhibits  a 
high  coherence  as  shown  in  Table  5.2.  The  observed  micropulsa¬ 
tion  has  a  high  degree  of  polarization.  Micropulsations  are 
generally  considered  to  be  left-hand  elliptically  polarized 
[Ref.  5].  The  observed  micropulsation  is  elliptically 
polarized. 

If  the  errors  in  phase  and  amplitude  that  were  shown  to 
exist  in  Section  IV  are  constant,  then  that  should  not  be  a 
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Figure  5.2.  Y  Coil  Magnetic  Field,  17  August  82,  1305-1310  Local. 

Amplitude  (nanoteslas  :  1  unit/in)  vs.  Time  (seconds  :  50  units/in) 


Amplitude  (nanoteslas  :  0.5  units/in)  vs.  Time  (seconds  :  50  units/in) 
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i  Table  5.2.  Observed  Micropulsation  Properties 


Frequency  (Hz):  0.058 


Peak-to-Peak  Amplitude  (nT):  Amplitude  =*yXJ  +  Y  +  Z2 


PSD  (dB)  ( Ref  lnT/*Hz) :  X-Coil:  -3dB 

Y-Coil:  -4dB 


z 

-Coil:  1. 5dB 

Coherence:  X-Y  Coils: 

0.99 

Y-Z  Coils: 

0.99 

Z-X  Coils: 

0.99 

Degree  of  Polarization: 

X-Y  Plane: 

0.99 

Y-Z  Plane: 

0.99 

Z-X  Plane: 

0.99 

Ellipticity: 

X-Y  Plane: 

0.45 

(left  hand) 

Y-Z  Plane: 

0.10 

(right  hand) 

Z-X  Plane: 

0.30 

(left  hand) 

serious  problem  since  we  are  dealing  with  a  single  frequency. 
Thus,  the  observed  micropulsation  exhibited  definite  structure 
decidedly  different  than  the  background. 


VI.  CONCLUSIONS 


A  type  PC-3  micropulsation  event  of  peak-to-peak  amplitude 
(0.37  nT)  and  period  (17.2  seconds)  was  identified  in  a  section 
of  data.  At  the  frequency  of  the  micropulsation,  the  coherence 
between  the  three  orthogonal  coil  sensors  was  about  0.99  for 
the  xy-pair,  0.99  for  the  xz-pair,  and  0.99  for  the  yz-pair. 

The  degree  of  polarization  was  0.99.  The  ellipticity  was  about 
0.30. 

The  origin  of  the  regular  "hump“-like  structure  in  the 
PSD 1 s  of  previous  coil  data  was  found  to  be  artificial  and  not 
representative  in  any  way  of  natural  phenomena  in  the 
geomagnetic  field.  Previous  conclusions  based  on  a  frequency 
domain  analysis  of  coil  data  characterized  by  a  regular  “hump1 
structure  (strong  or  weak)  in  the  PSD  plots  are  known  to  be 
unsubstantiated . 

Compared  to  a  PSD  frequency  domain  analysis  of  the  data 
set,  an  analysis  of  geomagnetic  field  fluctuations  in  the  time 
domain  easily  identifies  data  sets  containing  a  micropulsation 
event.  A  double  running  average  routine  developed  and  applied 
to  the  time  series  data  sets  was  found  to  be  helpful  in 
highlighting  the  micropulsation. 

Detailed  conclusions  about  the  geomagnetic  field  from  the 
coil  measurements  require  a  balance  of  analysis  in  the 
frequency  and  time  domain. 
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An  accurate  reproduction  of  a  time  series  plot  for  the 
geomagnetic  field  oscillations  will  require  modification  to  the 
existing  software. 

Comparisons  in  the  time  or  frequency  domain  between  any 
pair  of  coils  in  the  orthogonal  system  might  be  misleading  as  a 
result  of  cross  talk  between  electrical  inputs.  The  suspicion 
of  cross  talk  is  based  on  the  observed  pick-up  on  an  unloaded 
input  of  the  contents  of  an  adjacent  loaded  and  active  input. 

Further  micropulsation  studies  should  continue.  This 
initial  investigation  of  the  coherence,  polarization,  and 
ellipticity  features  of  micropulsation  events  shows  promise 
that  with  more  data,  one  might  suggest  novel  physical 
mechanisms  and  models  for  the  origin  of  micropulsations. 
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APPENDIX  A 


VOLTR  COMPUTER  PROGRAM 

Appendix  A  contains  the  VOLTR  computer  program  that 
provides  unsmoothed  time  series  voltage  plots.  VOLTR  uses  the 
subroutine  RD  to  strip  voltage  data  off  digital  tapes.  The 
integer  voltage  data  is  normalized  to  +  5V  and  plotted  using 
the  DRAWP  subroutine. 
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APPENDIX  B 

VOLTS  COMPUTER  PROGRAM 


Appendix  B  contains  the  VOLTS  computer  program  that 
provides  smoothed  time  series  voltage  plots.  VOLTS  uses  the 
subroutine  RD  to  strip  voltage  data  off  digital  tapes.  The 
integer  voltage  data  is  normalized  to  ♦  5V.  A  double  144  point 
running  average  is  applied  to  the  data.  The  DRAWP  subroutine 
plots  the  data. 
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APPENDIX  C 


LFVTC1  COMPUTER  PROGRAM 

Appendix  C  contains  the  LFVTC1  computer  program  that 
provides  smoothed  time  series  magnetic  field  plots.  LFVTC1 
uses  the  subroutine  RD  to  strip  voltage  data  off  digital  tapes. 
The  voltage  data  is  normalized  to  +  SV.  The  voltage  data  is 
fast  fourier  transformed  using  the  FOURT  subroutine.  The 
transfer  function  is  applied  and  forward  transform  taken.  The 
resulting  magnetic  field  data  undergoes  a  double  144  point 
running  average  and  is  plotted  using  the  DRAWP  subroutine. 
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